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ABSTRACT 



Satellite data first became available In the 1960*s 
making solar budget calculations more feasible. In this 
study, satellite measurements of reflected solar insolation 
from NIMBUS II were combined with simultaneous pyrheliometer 
readings of transmissivity over North America to provide the 
basis for computation of the atmospheric absorptivity. In 
the statistical parameterization of the absorptivity, such 
readily available data as water vapor mass and satellite re-- 
flectance were introduced as independent variables in the 
regression equations. Correlation coefficients between obser- 
ved regression-determined absorptivities and transmissivities 
were computed based upon an initial total of 235 data-samples , 
and the results interpreted. The conclusion was twofold: 

(1) A new suggested insolation- budget model is proposed for 
computation of absorptivity; (2) the desirability of deriv- 
ing a statistical parameterization for cloud-cover at grid 
points in terms of satellite-observed reflectances and window- 
channel temperatures is shown. 
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I. INTRODUCTION 



The solar energy cycle is of vast Importance to the 
meteorological community since addition of heat is a pri- 
mary source of energy in atmospheric circulations. There- 
fore, a proper accounting of absorptivity and transmissiv- 
ity within the atmosphere is vital to improved modeling 
and numerical forecast procedures. 

This study is designed to correlate satellite measure- 
ments of reflected solar Insolation from the Medium Resolu- 
tion Infrared Radiometer (HRIR) on NIMBUS II and simultan- 
eous pyrheliometer readings of transmissivity over North 
America. This research w as directed towards the parameter- 
ization of the atmospheric absorptivity and transmissivity 
in terms of available satellite albedo observations. In 
addition such other readily available data such as water 
vapor mass and cloud-cover may be introduced for perfect- 
ing the parameterization. 

The limits of the study were confined to the United 
States and Canada where data were readily available. 

Studies of the disposition of solar insolation prior to 
the acquisition of satellite data have been reported by 
Gabites (1950), London (1957) and Sellers ( 1965 ). Each of 
these researchers gave estimates of global albedo of close 
to 35,L. Vonder Haar ( 1968 ) using satellite data, has shown 
the global albedo to be close to 30 $. In view of this 
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reduction of the albedo of 5# from the previous estimates, 
one should look for revised transmissivity and absorptivity 
parameter izat ions which would account for the insolation 
balance . 

It would seem that successful parameterization of ab- 
sorptivity and transmissivity could be accomplished over 
data-rich areas such as existed over North America during 
the NIMBUS II period of this study. In the long-term, one 
would hope to describe these variables of the solar budget 
in terms of readily accessible satellite observations, and 
ultimately expand the parameterization to oceanic and other 
data-sparse areas. 
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II. DATA SOURCES AND PROBLEM OUTLINE 



A. GENERAL 

The data used in this study was from three basic sources 
and covered a five-day period from 15 July 1966 to 20 July 
1966. The surface data were compiled at the 47 pyrhelio- 
metric stations shown in Fig. 1. The locations were chosen 
to provide a broad coverage of solar transmission effects 
on the North American continent within the limits 25° north 
latitude to 59° north latitude. Several other types of 
data compilations (described below) were required to perform 
computations of solar reflectance, transmittance, and ab- 
sorptivity at pyrheliometric stations at the effective time 
of the NIMBUS II scan time. Thus all observational data 
had to be Interpolated to local solar times corresponding 
to the effective scan time of the NIMBUS II satellite at 
the pyrheliometric station. Solar time is computed every- 
where on the basis that 1200 solar time corresponds to the 
highest solar elevation of the day. 

B. NIMBUS II REFLECTANCE DATA 

NIMBUS II solar reflectance values were available for 
this study only in hard copy printout format. The NIMBUS 
II satellite had a nearly sun-synchronous polar orbit with 
the ascending node at 11:28 (equatorial solar time). The 
satellite scanned a swath perpendicular (Fig. 2) to its 
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Figure 1 - Mercator map-projection of data region showing 
pyrheliometer station locations used in this study. De- 
tailed latitude and longitude values for each station are 
located in Appendix Table A. 
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orbit, and extending between the horizons. Each scan spot 
is approximately a 30 -mile square at the subsatellite point 
and increases in size with increasing nadir angle of the 
scan spot relative to the satellite. 

The NIMBUS II Medium Resolution Infrared Radiometer 
(MRIR) was designed to scan in the five wavelength regions 
(listed below) with each region accompanied here by a 
brief description of the observational purpose of channel 
observation (e.g., NIMBUS II User’s Guide, 1966 ) . 

6.4 to 6.9 microns - This channel covers the 6.7 
micron water vapor absorption band. Its purpose is 
to provide information on water vapor distribution 
in the upper troposphere and, in conjunction with 
the other channels to provide data concerning rela- 
tive humidities at these altitudes. 

10 to 11 microns - Operating in an atmospheric 
’’window," this channel measures surface or near- 
surface temperatures over clear portions of the at- 
mosphere. It also provides cloud cover .and cloud 
height information (day and night). 

14 to 16 microns - This channel, centered about 
the strong absorption band of CO 2 at 15 microns, 
measures radiation which emanates primarily from the 
stratosphere. The information gained here is of pri- 
mary importance in following the seasonal stratos- 
pheric temperature changes. 
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LONGITUDE 




Figure 2 - Mercator map showing NIMBUS II orbital paths 
as indicated by the solid lines — . The limits of the 
swath coverage are indicated by the dashed lines — . 
Overlapped regions are indicated by the cross hatching//. 
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5 to 30 microns - This channel measures the emit- 
ted long wavelength Infrared flux and, In conjunction 
with the reflected solar radiation channel furnishes 
data on the heat budget of the planet . 

0.2 to 4.0 microns - This channel covers more than 
99 % of the solar spectrum and yields Information on 
the Intensity of reflected solar energy from the earth 
and Its atmosphere. 

These five channels are frequently termed channels 1, 

2, .... 5, respectively, of the MRIR experiment of NIMBUS 
II. Of the various channel-data available, only channel 5 
relating to solar reflectance from the earth and Its at- 
mosphere was utilized. 

It would perhaps have been a useful Indirect measure 
of cloud cover to have used the "window" channel-equivalent 
blackbody temperature data for the same observation times. 
However,- synoptic surface temperatures at the pyrheliometric 
stations would have been necessary to give a concise in- 
dication of the cloud amount and/or height, so that the 
channel 2 data were not analyzed here. 

The effective radiant reflectance R in the direction 
of the satellite from the reflecting surface filling the 
field of view of the radiometer and illuminated by unatten- 
uated solar radiation Is defined (p. 57. NIMBUS II Users’ 
Guide, 1966 ) by 
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H(I * , J' ,Z» ) = 



( 1 ) 



Sr_^ cos Z' 
m 



where I’.J’ are the coordinates of the scan spot, Z' Is the 
zenith angle of the sun and 



h = 



( 2 ) 



is the date-adjusted value of the solar constant S (correc- 
ted to 17 July 1966 ). The parameter N in (1) is the solar 

—2 —1 

intensity or radiance in watts m ster reflected back 
to the satellite, and is here assumed Independent of the 
viewing angles (nadir and azimuth) of the satellite rela- 
tive to the scan spot. The latter assumption is not ac- 
tually valid as has been shown by Raschke et al. ( 1973 ) 
who showed the distribution of the bidirectional reflec- 
tance over typical land spots to be dependent upon both 
angles. However, our scan-spot reflectances were computed 
by NASA according to Eq . 1. The hard-copy reflectance data 
obtained from NASA were in the form of averages composited 
to gridpolnts on a Mercator map projection of scale 1:10M 
(Fig. 3) and with a grid-mesh (I,J). This (I,J) grid-mesh 
may be formulated in terms of longitude A and latitude 0 
by means of the transformation: 

I - 1 = 0 . 8 ( 150 ° - A° ) (3) 
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Figure 3 - Mercator map-projection of data region used in 
analysis of data considered in this study. Equally spaced 
I- and J- intervals are indicated along the. bottom and 
left boundaries of the rectangular area. Corresponding 
longitude and latitudes are shown along the opposite boun- 
daries. 
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( 4 ) 



J - 1 = 0 . 8 In (I_±_§Hl£ 

7T COS rh 



COS <f> 



Some of the NIMBUS II data orbits used on 16 July 1966 are 
sketched in Fig. 2. 

Within any grid space centered at (I,J) the number of 
scan spots may vary from five to ten depending on the nadir 
angle to the scan spot, the actual earth-size of the grid- 
area (which is determined by the map factor m = sec <f> , dg = 
^map/ m ) * and the ext ent of swath overlap at the scan spot 
between successive orbits. 

The mean reflectance value obtained at each grid point 
(I,J) is therefore an average composited from all of the 
individual scan-spot reflectances R (I , J , I * , J ' ,Z * ) . Here 
I , ,J* are individual scan spots contained in the larger 
grid-mesh I,J of the map (Fig. 3) within the grid square. 
The averaged value R(I,J,Z) is the mean solar reflectance 
to space, computed after averaging over all the scan spots 
I',J' within I,J and is a function only of the grldpolnt 
identifiers (I,J) and of Z, the mean solar zenith angle at 
(I,J). R(I,J,Z) as thus determined is termed the direc- 
tional reflectance at the grid point, the direction being 
that of the mean solar zenith angle to the center of the 
grid square. Thus after the averaging reflectances to grid 
points has been accomplished, the resulting reflectance 
R (I , J,Z) is given by: 




( 5 ) 
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t 



where 

5 = solar constant = 2.0 langleys min -1 

cos Z = mean cosine of the solar zenity angle sub- 
tended by point I,J at station. 

cos Z = cos <f> cos 6 cos h + sin 4> sin 6 (6) 

<f> = station latitude 

6 = solar declination 

h = hour angle at effective time of satellite 

scan. 

N = mean reflected solar radiance from (I,J) to 

the satellite. 

(r m /r)= ratio of the mean distance to the sun and the 
distance to the sun at the required time. 

1 . Effective Tine Calculation 

The effective Greenwich Mean Time of the satellite 
reflectance observation at any station was computed using 
a linear interpolation scheme between the times along the 
swath over the station and orthogonal to successive satel- 
lite tracks, and the stations distances relative to the 
tracks as shown by Fig. 4. Thus Eq. (7) is used to inter- 
polate for effective time over the station as shown in 
Fig. 4. 

Effective scan-time = (^HTg) + (T-^ ) (7) 

If the station was in a region scanned only by one swath, 
the time T^ was read directly off the nearest track in 
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Figure 4. Interpolated. Scan-Time 
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Fig. 2. The effective time T was then converted to true 
solar time for the station by the following equation: 



S.T. 



T (GMT) 



(LONG. DEG) * 
15 (LONG. DEG) 



(HRS) 



( 8 ) 



Note the hour angle h depends upon the solar time of the 
scan as given by Eq. 9 and has the form: 



h = S.T. (hours) - 12 * 15° (9) 

2 . Station Reflect a nce Values 

The pyrhellometer station location (Fig. 1) falls 
within groups of four grid points shown In Fig. 5« There- 
fore, a linear Interpolation scheme was used to determine a 
value for the exact station location (I+4l, J+4J) located 
between the four adjacent grid points at which the reflect- 
ance values R(A), R(B), R(C), R(D) were available. The 
reflectance value at the location I+4l, J+4J (the station 
location) is given by means of the standard spatial inter- 
polation formula . 

R = (1-HJ) [(1-4I)R(A)+4I*R(B)] +A J [ (1-HI )R (D ) 

+HI*R(C)] (10) 

Here R(A) Is the reflectance at point A, R(B) Is that at 

B, etc., as Indicated by the grid array In Fig. 5« 

C. SURFACE TRANSMITTANCE OF SOLAR INSOLATION 

Surface transmittance of solar radiation was based on 
pyrhellometer readings recorded at the stations In Fig. 1. 
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Reflectance-Interpolation to Stations. 
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These records provided on magnetic tape are listings of the 
hourly accumulations of solar Insolation observed at the 
end of each true solar hour. The taped record was kindly 
furnished with the technical assistance of the Commander, 
Naval Weather Service. The values listed were pyrhelio- 
meter readings of the total hemispheric insolation, with 
the pyrhellometer calibrated to conform to a solar constant 
of 2.0 langleys min“l. 

The locations of pyrheliometric stations were converted 
from latitude 0 and longitude X to the (I,J) coordinate 
system in accordance with Eqs. 3 and 4. The mean station 
transmittances were determined at the effective time of the 
satellite scan. This was done to give a transmittance 
coincident with the "simultaneous" satellite mean reflect- 
ance observations over this station. The actual station 
insolation was considered to be time-centered at 0930 solar 
hours when it was reported for the solar hour ending at 
1000, etc. Thus, it was necessary to interpolate an in- 
solation value at the station coincident with the effective 
time of the satellite-scan. This was accomplished by using 
a cubic Lagranglan time- interpolation scheme described 
below. 

The four successive values of the hourly insolation that 
bracketed the effective scan time were used in the following 
equation: 
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Here 



Qg - 



(T-r2H(T-n)*(T-r4UQi 
(T1-T2 )* (T1-T3 ) * (Tl-1’4 ) + 



(T-Tl ) * (T-T3 ) » (T-I4 ) *92 
(T2-T1 ) * (T2-T3 ) * (T2-T4 ) 



+ 



(T-Tl)«(T-T2)»(TrT4)»33 . (T-Tl )* (T-T2 ) * (T-T 
(T3-T1 ) * (T3-T2 ) * (T3-T4) + (t4-T1)* (T4-T2 )* (T 




( 11 ) 



Q Q = station insolation interpolated to the 
scan time. 

T = effective time of satellite-scan. 

T1,T2,T3,T4 = successive times of pyrheliometric 

readings at half-hour solar clock times. 

Q1.Q2,Q3.Q4 = station insolations at the four suc- 
cessive times T1,'T2,T3.T4. 

The Langrangian fit was used to reduce the effect of 
transistory clouds over a short time frame. 

D. PRECIPITABLE WATER VAPOR 

Precipitable water vapor was computed from radiosonde 
soundings nearest to the pyrheliometric stations. The 
soundings were taken at 0000GMT and 1200GMT starting at 
OOOOGMT 16 July 1966, and continuing through 0000GHT 20 
July 1966 . 16 levels were employed using the surface as 

the bottom level and the 15 standard levels at whole mul- 
tiples of 50 mb above the surface. Temperatures and rela- 
tive humidities (RH) at these levels were extracted from 
the soundings and used to compute water vapor mixing ratios 
at the indicated levels in the vertical. The computation 



25 



of q was facilitated, with RH known, by use of the Goff- 
Gratch (cf., List, 1958, P» 3 50) formulation of the satura- 
tion vapor pressure 



log 1() e w = -7. 90298 (T s /T-l)+5. 02809 log 10 (T s /T)- 
1.3816 X i0-7(10ll-344(l-T/T s ).i) + 



8.1328 * 10-3 * ( 10 -3.49149(T s /T-1L 1) + 



los 10 e ws 



( 12 ) 



where 



e w = saturation vapor pressure over a plane surface of 
pure ordinary liquid water. 

T = absolute temperature (°K). 

T g = steam point temperature (373*l6°K). 
e ws = sa ^uration pressure of pure liquid water at steam 
point temperature (1 standard atmosphere = 
1013.246 mb). 

Following the use of Eq. (12), the value of q is given by: 



Then, precipitable water vapor is found from the mixing 
ratio integration 



q = 621.97 = 621 . 97 (RH)e 

p-e — 



p-(RH)e w 



w 



(13) 




( 14 ) 
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For the times on 16 July 1966 through 19 July 1966 a linear 
interpolation scheme between the 1200 and 0000GMT soundings 
was used to establish the preclpltable water vapor at the 
effective time of satellite scan at each station. 

More specific examples of the calculations performed 
in order to relate the satellite reflectance, and simul- 
taneous stat ion-transmittance through a physical relation- 
ship are detailed in Section III. 
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III. DATA PROCESSING 



A. THE SOLAR-INSOLATION MODEL 

In Section II, a general description was given of the 
data-collection problem at each pyrheliometric station at 
the effective time of the satellite scan. The immediate 
objective of the data-analysis corresponding to each station- 
observation itfas to solve the solar- insolation budget equa- 
tion for the atmospheric absorptivity. 

FIGURE 6 




QrI Top of 
atmosphere 



Qa 



^gI a G"-of 

E5HTET STJHFACE 

Figure 6. The solar-insolation model 
over a pyrheliometric station. See 
Eq . 15 and following for explanation 
of symbols. 



Fig. 6 shows an insolation model after Fritz et al. (1964) 
and Hanson et al. (196?) which Illustrates the terms involved 
in the solar-insolation budget at a satellite scan-time for 
a given pyrheliometric station. These terms ares 

Q s = S(r m /r) 2 cos Z (15) 
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Q s = Incident extraterrestial insolation based upon a 

solar constant S of 2.00 langleys min”^ (Qg is the 
denominator of Eq . 5)* 

Qq = pyrheliometric value of the insolation at the sta- 
tion interpolated to the true solar time (Eq. 8) 
of the satellite scan, according to Eq. 11. 

Qr = RQs * which is the reflected insolation at the sta- 
tion computed using the space-interpolated reflect- 
ance R(I+zlI, J+dJ) of Eq. 10. 

= absorbed insolation in the atmospheric column above 
the station. 

a G = surface albedo of the earth in the vicinity of the 
station. 

The midsummer values of c* G were taken from a North 
American atlas reproduced here as Fig. 7 (prepared by Rung 
et al., 1964). The term l-a G is the fractional absorption 
of the incident insolation ^ , at earth. At each of the 4? 
pyrheliomet er stations used, the atlas value of a G has been 
listed next to the station coordinates (Appendix A). 

The insolation-budget equation then according to Hanson 
et al. (1967) becomes: 

Qg = Qpj + Q a + Qq ( 1 ~ ct q ) (16) 

If each term of (16) is divided by Q s , the equation may be 
rewritten as 
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Figure 7 



Summer surface albedo map for North America 
(after Rung et al., 1964) 
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q A 1 ” “ q G a G^ 



(17) 



where : 

q A = atmospheric absorptivity. 

<1r = Qr /^3 Is the computed reflectance using the NIMBUS 
II data-maps. 

qQ = atmospheric transmissivity Qq/Q 3 as inferred from 
the pyrheliometer reading. 

The parameter a G was taken from the summer mean reflect- 
ance maps of Kung et al . (1964) for each station, and assum- 
ed constant during the NIMBUS II period 15-19 July 1966, al- 
though a„ varied from station to station in accordance with 
G 

varying terrain and vegetation characteristics. The use of 
the a G in connection with Sq . 1? was considered valid by 
Hanson et al. (1967) both in clear and partly cloudy sky 
conditions, and this same approximation is made here in 
reference to the computation of q^ in conditions of mixed 
cloudiness. An attempt is made in the subsequent regression 
analysis to isolate that part of q^ which depends upon 
cloudiness. Statistical specification of q a similar to that 
of Hanson will be performed here using the simultaneously 
observed and/or calculated values of q^, q G , and ot during 
the period 15-19 July 1966, with q R measured by NIMBUS II 
satellite . 

B. DATA TABULATIONS 

For present tabulation purposes, -the 47 insolation sta- 
tions were assigned arbitrary station numbers 1, 2, ..., 47 
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TABLE I. TYPICAL SET OF INSOLATIONAL 
AND RELATED DATA AT PYRHELIOriETER STATIONS AT 
THE EWffiTTO TIKE OF THE NIMBUS II SCAN ON 17 JULY 1966 



Station 

No. 


Extra- 

terrest 

Qs » 

ly nin~ x 


Re- 

flected 

q r 


Pyrhel- 

lometer 

Qq 


Absorbed 

(Air) 

Q. A 


Preclp 

Water 

Vapor 

U,gm cm” 2 


Zenith 
Path 
SEC Z 




1 


1.7790 


0.6398 


0.8319 


0.4404 


4.8232 


1.0883 




2 


1.8857 


0.2711 


1.3564 


0.4210 


4.5175 


1.0267 




3 


1.9189 


0.3262 


1.6974 


0.1500 


4.3973 


1.0089 




4 


1.7939 


0.3243 


1.4725 


0.2328 


4.4853 


1.0792 




5 


1.8884 


0.7836 


0 . 6680 


0.5236 


4.3882 


1.0252 




6 


1.8381 


0.6659 


0.6116 


0 . 6461 


4.3088 


1.0533 




7 


1.7782 


0.3921 


1.0544 


0.5004 


3. 0089 


1.0887 




8 


1.7949 


0.3647 


1.1184 


0.5130 


2.6633 


1.0786 




9 


1.6741 


0.2470 


1.2202 


0.3778 


1.8045 


1.1565 




10 


1.7562 


0.3273 


1.1706 


0.4573 


2.2867 


1.1024 




11 


1.8971 


O.4945 


1.1978 


0.4324 


2.4511 


1.0205 




12 


1.8742 


0.4056 


1.6545 


0.1782 


2.1694 


1.0330 




13 


1.8047 


0.3676 


1.3676 


0.3294 


1.3708 


1.0727 




14 


1.8233 


0.2287 


1.3959 


0.4080 


2.1878 


1.0618 




15 


1.8660 


0.3730 


1.2186 


0.5424 


3.2688 


1.0376 




16 


1.8013 


0.2279 


1.3395 


0.4348 


1.3032 


1.0748 


P 


17 


I.6696 


0.5659 


-0.0474 


1.1439 


3.0095 


1.1596 




18 


1.7226 


0.3008 


1.1724 


0.4487 


1.5689 


1.1238 




19 


1.7201 


0.2449 


1.3934 


0.3048 


1.6416 


1.1255 


A 


20 


1.6415 


0.5903 


1.3108 


-0.0499 


1.6987 


1.1794 




21 


1.7945 


0.3556 


1.4354 


0.1901 


1.6033 


1.0789 




22 


1.8787 


1.0111 


0.2224 


0 .6764 


4.5055 


1.0305 




23 


1.7861 


O.38O6 


1.2782 


0.3191 


2.4676 


1.0839 




24 


1.7012 


0.2461 


1.0973 


0.5225 


1.4247 


1.1380 




25 


1.7691 


0.3044 


1.2806 


0.3762 


1.5^37 


1.0943 




26 


1.7277 


0.3217 


1.3567 


0.2528 


3.1032 


1.1206 




27 


1.5645 


0.5559 


0.3864 


0.6763 


2.2638 


1.2375 


A 


28 


1.4358 


0.7002 


1.0219 


-0.1330 


2.1568 


1.3484 




29 


1.6297 


0.5130 


1.1258 


0.1598 


1.6994 


1.1879 




30 


1.5447 


0 . 4654 


1.0314 


0.2027 


2.2699 


1.2533 




31 


1.6977 


0.2495 


1.2607 


0.3766 


1.8045 


1.1404 




32 


1-5499 


0.4177 


1.0461 


0.2535 


2.1985 


1.2491 




33 


1.7278 


0.2937 


1.2433 


0.3648 


1.4729 


1.1204 




34 


1.7287 


0.2173 


1.2437 


0 .4666 


1.3883 


1.1199 




35 


1.6254 


0.5605 


0.4457 


0 .6861 


3.0243 


1.1910 




36 


1.7199 


0.2183 


1.1793 


0.4992 


1.9092 


1.1256 




37 


1.6432 


0.2796 


1.2288 


0.3315 


2.1046 


1.1781 




38 


1.6118 


0.6234 


0.7737 


0.3308 


2.8319 


1.2012 




39 


1.6971 


0.3342 


1.1203 


0.4218 


2.1268 


1.1408 


P 


40 


1.6605 


0.1193 


1.6500 


0.1553 


1.4620 


1.1659 


P 


41 


1.6035 


0.3091 


1.6500 


-0.1081 


1.5630 


1.2074 




42 


1.6532 


0.5877 


1.7551 


0.4311 


1.5231 


1.17H 




43 


1.6297 


0.2803 


1.1822 


0.3445 


1.6260 


1.1879 




44 


1.7255 


0.2697 


1.2719 


0.3620 


1.4731 


1.1220 




45 


1.7255 


0.2697 


1.2948 


0.3423 


1.4731 


1.1220 




46 


1.6361 


0.4925 


0.9245 


0.3671 


1.6994 


1.1833 




47 


1.6160 


0.6268 


0 . 6268 


0.4459 


3.1851 


1.1966 
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(see Appendix A, for station names, code and geographical 
coordinates). Thus, for each station of Fig. 1, the values 
Qg, Q q , were computed using Eqs. 15. 11. 10; and 

was computed as the residual from Eq . l 6 . A typical set 
of insolation values obtained for each of the 47 stations 
on 17 July 1966 is listed in Table I. 

The primary objective here at the outset was to obtain 
a statistical regression for the absorptivity q^ against 
an effective measure of the water vapor absorbing mass in 
the air column above the pyrhellometer station. For this 
purpose simultaneous values of preclpltable water vapor U 
and of sec Z were also listed in Table I. 

The complete set of raw data to be used for this re- 
gression analysis is tabulated in Appendix B, There were 
in all 47 X 5 = 235 simultaneous values of the variables 
Qg. Qr» Qq* Q a . as these have been defined above in 
Eq. 15 and in the lines just following- Sq. 15. Of the 188 
reports listing the dependent sample-data taken from days 
2 through 5 (July 15. .... 19th, inclusive) in Appendix B, 
it was found that eight cases had missing pyrhellometer 
values at the solar hour adjacent to the effective time of 
the NIHBUS scan, and thus led to unreasonable transmissivit- 
ies. These pyrheliometer-assoclated inconsistencies in the 
data field were denoted by the code symbol P. 

On the other hand, for days two through five, there also 
occurred six cases of negative after Eq . 16 was applied, 
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which were also physically impossible. All cases for which 
Q A were negative were coded, by the symbol A (negative ab- 
sorptivity). Hence the remaining sample of 174 data sets 
from days two through five were used in the further screening- 
analyses to be described in Section III (C). For an illus- 
tration of the initial screening involving data-deletions 
of code P and A, reference should be made to Table I in 
the column listing the pyrheliometer station number. Sec- 
ondary data-screens which lead to further data-exclusions 
based upon statistical tests, are utilized for additional 
deletions of data cases based upon improbable extreme values 
of 



Q A ’ 9 R 



These tests are essentially based upon the requirements: 



q A = 1 



q G (1-Cfc) - q R 



<q A (MAX) 
>q A (KIN) 



( 18 ) 



Here the values of q^(I>iIN) and q^(MAX) are indicated by 
empirical evidence to be discussed in Section III (C). 



C. APPLICATION OF A DATA -LI MIT SCREEN 
1 . General Considerations 

Hanson et al . (1967) developed a best-fit regression 
for q A using data based upon TIROS reflectances, and an em- 
pirically-determined measure Ify of the water vapor absorber 
mass. Hanson found the most effective definition of h^, which 
led to the highest linear correlation of q^ with M^, to be: 
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(19) 



M w = (U sec 2)* In (U sec Z) 

The specific form of Hanson's regression was: 

q A = 0.096 + 0.045 M h . (20) 

which in turn was optimal in yielding the highest value of 
the correlation coefficient R between q^ and My 

R (q. . Ii ) = Maximum 
A V/ 

Here My of (19) was one of many possible functions of U sec 
Z having the empirical form of an absorptivity (see Fig. 8). 

Application of the HIKED 02R statistical program (on 
file in the Program Library of the W. R. Church Computer 
Center of the Naval Postgraduate School) to the dependent 
sample of 174 data-cases from days two through five, was 
made in order to generate the least-squares best-f it between 

(a) q A and My 

(b) q G and My 

(c) q R and My 

The most useful immediate results obtained here were 
listings of sample means and standard deviations (S.D.) of 
each variable under cons ideration, (as well as listings 
of the 174 values) of q A , q^ and q R (see Table II). 
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Figure 8 . Absorptivity of solar radiation In the atmos- 
phere as a function of optical pathlength. In U Sec Z. 
Curve (1) is the absorptivity in a clear -atmosphere by 
H 2 O, CO 2 . and O 2 t from Yamamoto ( 1962 ). Curve (2) is 
the absorptivity in the clear atmosphere by H 2 O only, 
from Houghton (195*0 • Curve (3)- is the absorptivity in 
the atmosphere over the United States with mixed cloud- 
iness, based on the Hanson's 19&7 study. Vertical bars 
indicate the most probable error at the abscissa value. 
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TABLE II. 


Statistics 


Pertaining to q A , q G , q p 




Mean 


S.D. 


Range 


q A 


0.2290 


0.0943 


0.0026 to 0.6192 




0.6308 


0.1964 


0.0747 to 0.9594 




0.2388 


0.1127 


0.0916 to 0.6406 = qg(i-iAX) 



2 . Conslderat Ion of Data-Kxcluslons on q^ 

The listed range of the planetary albedo q^ in 
Appendix B is within the range of albedo values quoted in 
the literature (List, 1958) for various surfaces: cloud- 

tops, soil-types, etc. Hence no limiting range on the values 
of q^ found here were imposed on any of the data cases of 
this study, i.e., no data-case exclusions were made on the 
basis of q R . This does not mean that q^ does not have a 
measurement error, but that its effect on upsetting the 
balance of Eq. 17 is not as great as the more unrealistic 
ranges of values associated with q A and q G . 

3. Data-Exclusions Based on q. 



Because q A was determined as a residual involving 
the two measured (but error-subject) quantities q~ and q G , 
it was felt that some reasonable bounds would have to be 
applied to the range of q^. Thus q^ (MIN ) and q A (MAX) values 
were established in accordance with extreme-probability 
limits on both wings of a Gaussian distribution centered on 
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q^(MEAN). From the statistics of Table II, It follows that 
the definition of: 

q A (MIN) = 0.2290 - 2.3267 (0.0943) = 0.0096 (21) 

q A (MAX) = 0.2290 + 2.326 7 (0.0943) = 0.4484 (22) 

The value of q A (I*IIN) was relaxed to: 

q A (MIN) = 0.07 (23) 

[of- Fig. 8, after Hanson et al. (1967)J* With the limits 
established by Eqs . 22 and 23 , there were six case-deletions 
due to small absorptivity-calculations 

q A < 0.07 

and 4 case deletions due to unacceptably large q, 

A 

q A > 0.4484 

These ten deletions are coded A MIN, A MAX In Appendix 
Table B. The codes KIN and MAX Indicate that statistical- 
empirical considerations were Involved In these q A secondary- 
deletions . 

4. Data-Excluslons based on qQ 

Here two statistical relationships analogous to Eqs. 
21, 22 for considering probable limits on q^ were set up: 



qQ (MIN ) = 0.6308 - 2.576 


(0.1964) = 0.1249 


(24) 


q Q (MAX) = 0.6308 + 2.576 


(0.1964) = 1.1367 


(25) 
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The factors 2.576 in £qs. 24, 25 correspond to Halting data- 
cases lying In a 0.5^ rather than in a 1# Gaussian probabil- 
ity distribution wing. Here we are being slightly more 
cautious about rejecting extreme values of q G than in the 
corresponding consideration of q^. Finally Eq. 25 has been 
altered from an extreme transmissivity of ll4;S to 93^ 

(l.e., 0.93). This follows in accordance with the minimum 
q^ = 0.07 of Eqs. 23 and 18 so that: 

q G (MAX) =0.93 (26) 

the upper limit exclusion Eq . 26 accounted for no data- 
case rejections, while the lower limit of Eq. 24 accounted 
for four data-case rejections. These cases are listed in 
Appendix Table B by the symbolism T MI W . The revised data- 
sample size for regression purposes now comprises 160 cases. 

Actually the purpose of reducing the number of extreme 
data cases was to eliminate as much noise as possible from 
the regression analysis so that the physical-statistical 
relationships would stand out more clearly. It will be 
seen in Section IV that the data-screenlng process for 
elimination of extreme-cases did not improve the regression 
analysis markedly, after the initial reduction to the sample 
of 174 cases (based upon obviously missing and erroneous 
data) was achieved. 
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IV. THE REGRESSION RESULTS 



A. THE REGRESSION ANALYSIS FOR q A 

Initially a least-squares fit analogous to that of 
Hanson (Eqs. 19. 20) was sought, using both the full sample- 
set of 174 cases, and the smaller sample-set of 160 cases 
after extreme-value limits had been applied to screen 
statistically for possible unreliable data. The regression 
method used was that of 3IMED 02R (Dixon, 1970). Seven 
variables were listed on an IBM card for each data-case. 
Following are the variables listed for each case in Appendix 
Table B, namely: 



Qs * Qg * Qq • » U , cos Z , Of q 

By means of transgeneration (an option of the BIMED series 
of programs), the fractions q A = Q^/Qg* Qq * were formed, 
as well as the water vapor mass (Eq. 19) used by Hanson 
et al. (1967). 

Two types of regression tests were run for both the raw- 
data and screened-data sample sets. These tests involved: 

(1) q A on M ^ alone 

(2) q^ on both M^- and q^ 

The reason for performing the regression-test (2) Is that 
it was observed that an above (below) normal q A tended to 
occur with an above (below) normal q^. In fact the correla- 
tion coefficient between q A and q^ was R(q^, q^) = 0.335* 
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Wore specifically, an Increase In q^ Is, In general, as- 
sociated with an Increase of cloud cover, so that the In- 
clusion of q R as a predictor for q A should Indicate at 
least partially an effect of the cloud cover on q^. More- 
over with the NOAA satellite scanning radiometer readout 
available at Fleet Weather Centrals, values of q R ought to 
be as readily accessible as the radiosonde data used for 
the M. computation. 

The two regression tests on the raw and statistically- 
screened data sets led to the following regression equations 
and associated correlation coefficients R. 

Test 1, sample size n = 174 (raw data test): 



q A = 0.1983 + 0.0188 

R = 0.1807 



( 2 ?) 



Test 1, sample size n = 160 (screened data): 
q A = 0.1978 + 0.0189 M w 
R = 0.2206 



(28) 



Test 2, sample size n = 174 (raw data): 
q A = 0.1649 + 0.1880 q R + 0.0118 M w 
R = 0.2802 



(29) 



Test 2, sample size n = 160 (screened data): 



q A = 0.1605 + 0.2162 q R + 0.0111 n w 
R = 0.3567 



(30) 



4l 



2 

Since the correlation coefficient squared (H ) Is the 
measure of the explained fraction of the variance Intro- 
duced by the regression technique. It Is clear that q^ Is 
not well specified by either ii^ and/or by q^ in spite of 
evidence to the contrary by Hanson et al. ( 1967 ). The 
most significant correlation for q^ of the set of Eqs. 

27 , 28, 29 , 30 , Is that Involving both q.^ and as pre- 
dlctands with statistically screened data. The mean sta- 
tistics which relate to the four tests Just cited are 
summarized In Table III. The mean standard error of es- 
timate for the best case (line 4) In Table III was S.E. = 

0.0715. 

The relatively large standard error which remains after 
the regression tests (Table III), seems to be due to the 
computational formula for q^ of (17). 

«A = 1 ' q H - (1 - a G )q G <17) 

This equation has been used with partly cloudy to cloudy 
skies, when q^ has already been reduced by the cloudiness 
associated with the large reflectance q^. However, In 
these cloudy cases, the multiplicative reduction of the q^ 
term of (17) by the factor (1-a^) gives a positive bias 
to the computation of q^ which shows up In the predominance 
of large positive residuals • remaining after application of 
the regressions, Eqs. 27, 28, 29, 30. 
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Table III 



Regression Statistics for 
Tests of (1) q A on M w ; (2) q A on , q R 

Standard 





Test 




Sample 

S},ze 


Mean Deviation 
QA 0 (Qa) 


Corr . 
Coeff. 


% ExpL 
^krlance 


Std. Err 
of Est . 


(1) 


q A 


on 


“w 


1?4 


0.2290 


0.0943 


0.1807 


3.27 


0.0931 


(1) 


q A 


on 


M w 


160 


0.2275 


0.0760 


0.2206 


4.87 


0. 0744 


(2) 


q A 


on 


%* q R 


1?4 


0.2290 


0.0943 


0.2802 


7.85 


0.0911 


(2) 


q A 


on 


M tf * q R 


160 


0.2275 


0.0760 


0.3567 


12.72 


0.0715 
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B. A REGRESSION ANALYSIS FOR q u 

The relatively weak statistical specification of was 
due in part to improper accounting for earth-surface absorp- 
tion term in Eq. 17 when clouds were present. However q G 
is a measured parameter, which is related to q^ through 

Clear case: q G = (31) 

1_a G 

Overcast case: q G = l-q^-q„ (32) 

Hence, since the computational error in determining q^ tends 

to obscure the statistical regression of q A on My and q. . 

Eqs. 31. 32 suggest that such a statistical relationship 

might be contained in a regression of the transmissivity 

q G on My and q R in view of the functional relationships (3D. 

(32) between q~ and q.. There seems to be a greater likeli- 
u A 

hood of a significant regression for q G than for q^ in view 
of the fact that q is directly measured. The assumption of 
a significant regression of the type just described is the 
hypothesis of subsection IV (B). 

In the statistical tests to be conducted here, q is 
established as the predlctand in a set of regression tests 
for 

(1) q G on My alone, 

(2) q G on My and q^ . 

These tests then led to the following regression results in 
Eqs. 33 » 3^» 35. 36 which are the analogs of Eqs. 27. 28, 29 , 
30, respectively. The associated regression statistics are 
listed in Table IV. 
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Table IV 



Regression Statistics for 



Tests of (1) q G on 





Test 


Sample 

Size 


i v Iean 

q G 


( 1 ) 


sg 


on 


M W 


174 


0.6308 


( 1 ) 


q G 


on 


M W 


160 


0.6429 


( 2 ) 


q G 


on 


M w .q R 


174 


0.6308 


( 2 ) 


q G 


on 


M W q R 


160 


0.6429 



( 2 ) q G on M w and q R 



Std . 


Corr . 
Coeff . 


% Sxpl. 
Variance 


Std. Err. 
of Est. 


0.1964 


0.3124 


9.77 


0.1871 


0.1752 


0.3369 


11.35 


0.1654 


0.1964 


0.6236 


67.83 


0.1120 


0.1752 


0.8679 


75.32 


0.0676 
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Test 1, sample size n = 174 (raw data): 



q G = 0.74141 - 0.0 677M W (33) 

R = 0.3124 

Test 1, sample size n = 160 (screened data): 

q G = 0.74743 - 0.066 5My (34) 

R = 0.3369 

Test 2, sample size n = 174 (raw data): 

q G = 0.9887 - 1.392?q H - 0.0155^.; (35) 

R = 0.8236 

Test 2, sample size n = 160 (screened data): 

q G = 0.99252 - 1.42078q R - 0.01509^ (36) 

R = 0.8679 

Note the relatively high correlation coefficients in 
the Test (2) case in the specification of q G in terms of 
both M w and q„. The results seem to be slightly better in 
the case of filtered data (sample size of 160) , where the 
standard error of estimate of q^ is O.O 876 compared to 
0.1120 in the case of the raw data. 

In Table IV, the lower values of the correlation of q G 
on alone suggests that the q G estimations developed from 

(33). (34) will not be as useful as those derived from (35). 
( 36 ). This will be further explained in Subsection (C ) . 

C. TESTS ON INDEPENDENT DATA-SAMPLES 

In Sections IV (A) and (b) , four sets of regression 
equations for both q,, and q were statistically deduced 
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based upon the dependent -data samples of 174 and loO cases, 
respectively (termed the raw-data and statistically-screened 
samples, respectively). 

The regression equations so derived were of the form: 



and represent approximations to the observed or measured 
q^, q G* ^ qs * 27 and 28, the special case 

A x = 0 

was enforced, while in Eqs. 33. 34, we required 
B 1 = 0 

Eqs. 37. 38 with A^ = 0, = 0, respectively were termed 

the Test 1 regressions, while A^ ^ 0, ^ 0, have been 

termed the Test 2 regressions. 

The generating formulas of Eqs. 37 and 38 may be written 
in the matrix-product notation: 



^A = A 0 + A l q R + VV 



(37) 




(38) 




(39) 




(40) 



Here A Q , A lf a 2 are listed in Eqs. 37, 38, 39, 40 and B , 
, B 2 are listed in (33 36 ). 



• • • • 
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Table V 





Independent Sample 


Regression 


Tests 






(a) Raw-Data Sample 


, n=44 










Predic- 

tand 


Means 


Standard 
Dev. a(q) 


Standard 

Error 


Correlat Ion 
Coefficient 




q A =0.21?9 


0.0882 










q A 












t A =0.228l 


0.0263 


0.0860 


R= 0 . 2652 , 


by 


Eq . 27 




^= 0.2268 

q G =0.65?5 


0.0290 

0.1718 


0.0871 


R=0.2155, 


by 


Eq . 29 


q G 












f G = 0 . 63?2 


0.0805 


0.1688 


R=0 .2387 . 


by 


Eq. 33 




^ G = 0 . 646 5 


0.1544 


0.1073 


R=0 . 7868 , 


by 


Eq. 35* 


(b) 


Screened Data 


, n=4l 












q A = 0.2124 
f A =0.2253 


0.0688 

0.0198 










q A 


0.0688 


R=o.055i, 


by 


Eq. 28 




q A =0.2240 

q G =0.6745 


0.0246 

0.1467 


0.0688 


R=o.0975, 


by 


Eq. 30 


q G 












q“ G =0 . 6508 


0.0696 


0.1467 


R=0.0865, 


by 


Eq . 3^ 




^= 0.6588 


0.1383 


0.0882 


R=0 . 8048 , 


by 


Eq. 36* 
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For the independent test, the simultaneous samples 
(q A . q^) and (q G , *q G ) from day 1 were used. For the test 
on the raw-data independent sample of day one, the four 
specific coefficient sets of Aq, A^, Ag and 5 q, B^, of 
Eqs. 27, 29. 33* 35 were used. For the test on the screened 
data independent sample of day one, the four coefficient 
sets of Eqs. 28, 30, 34, 36, respectively were used. The 
regression results for the eight combinations of the one- 
predictor and two-predictor formulas for q, and q„ are shown 
in Table V for both the raw and screened independent sample 
tests. 

Table V shows that the use of screened data in place of 
raw data caused a decrease in the standard deviation of the 
predictand by eliminating some of the extreme balues. More- 
over the correlation coefficients: 



R(q A . “5^ • R(c k;* ^ 

decreased from the raw to the screened cases when the listed 
values of the correlation coefficients of the initial spec- 
ification Eqs. 27, 28, 29, 30, 33* 34 were small (in the 
range .2 to .35)* However the two key results which stand 
out in Table V are those corresponding to the use of specify 
ing Eqs. 35* 36 (which gave R = 0.7868 and 0.8048, respec- 
tively), denoted by an asterisk superscript in Table V. 

These two key results indicate that the transmissivity- 
estimator 'q'Q expressed by either Eq. 35 or 36 gives a sig- 
nificantly good specification for the pyrheliometric obser- 
vation of q G . 
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The interpretation of the correlation coefficients R 
in the two cases noted in Table V may be seen from the 
standard statistical formula: 

^2 _ Mean square explained by Regression 
Mean square of predictand 




Thus by use of Eq . 35. the standard error is 0.10?3, 
whereas the original standard deviation is 0.1?l8. Like- 
wise by use of Eq . J6 the standard error is 0.0882 compared 
to the original standard deviation of 0.146?, as summarized 
by Table V. 

Actually, there has been no appreciable significant 
change in the fractional part of explained variance by (35) 
versus (36). Moreover the process of providing screening 
limits is very laborious if it provides no significant im- 
provement in the fractional explained variance (right side 
of Eq . 4l ) . 

D. ERROR CONSIDERATIONS 

As has already been noted, there is a substantial stand 
ard deviation in q^ (0.1964 in the case of raw data, 0.1752 
in the case of screened). A substantial part of this error 
is due to Initial calibration of the Eppley pyrheliometers , 
ranging from 3*5 to 10.5/k» Hanson et al. (19?2). While 
such errors may be partially accounted for and correction 
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factors applied, the nature of the correction process was 
not known until the present study was completed. It Is 
likely that a sizeable fraction of the standard error of 
estimate of q G listed In the last line of Table IV is due 
to the pyrheliometer calibration error. 

Another minor source of instrumental error is that due 
to the necessity of interpolation of q G In time bet ween 
tabulated hourly solar transmission increments. The non- 
availability of detailed digitized records of the pyrhelio- 
meter and satellite reflectance data both on a sensitive 
time scale was another error-source in this study since 
short-term cloud variations of a period less than a half- 
hour would have been related to the short term variabilities 
of q G and of q^. Note that in this connection, our time- 
smoothed records of q G and q^ were highly correlated, there 
being a correlation coefficient R (q G , q- t ) = -0.821. 

As mentioned earlier in connection with the effect of 
cloud-cover on the computation of q G (Eqs. 31. 32 ). two 
different formulas apply in the clear and overcast cases 
for the solar-insolation budget as accounted for at the 
pyrheliometer. Failure to explicitly account for the amount 
of cloud-cover has biased the computation of q^ in the sense 
of increasing the value of the latter by up to 15/2 (that 
is, by the factor CL/(l-a G ) in Eq . 31. where CL is the cloud- 
cover). A discussion for the systematic inclusion of the 
cloud cover parameter for improved specification of q^ and 
q G is given in Section V. 
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V. RECOMMENDATIONS AND CONCLUSIONS 



As deduced from this study, the insolational budgetary 
model of Eqs. 1? and 31 gave q^ values biased towards high 
values. This was due to the failure to account for the 
effects of cloud cover properly, as shown for the overcast 
case in Eq. 32. The combination of Eqs. 31. 32 show that 
for some intermediate cloud-cover state CL, the properly 
weighted form of these two equations leads to 

q A = 1 " q R " q G t 1 " a G ( 1 - CL )] = 0 (**2) 

In Eq. 42, the term involving is the fractional insola- 

tion which the cloud-openings (1-CL) permit to escape to 
space. 

In this study, CL was not known and was assumed zero, 
as was also done by Hanson et al . ( 1967 ) and Fritz et al. 
(1964). However, recent studies by Downey et al . (1972) 

showed that q^ is highly correlated with CL, as indicated 
by the following summary of his results (based on NIMBUS III 
observations ) . 

Correlation coefficients R(q^, CL) by periods 1969-70. 

May 16-31 Aug. 1-15 Oct. 3-17 21 Jan-3 Feb. 

,, . ■ , I . 970 

5764 0.78 0.78 0.85 

With results of this type, combined with satellite 
measurements of window-channel temperatures as well, seasonal 
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regressions of some physically realistic form should be 
fitted to the observed cloud-coverage. The objective is 
to obtain a highly specifying representation for CL of form 

CL = F (q R , T bb ) (43) 

where F is the desired statistically-generated function of 
the indicated variables. 

A seasonal variation may be built into the function F 
by making F dependent upon sec Z, where Z = zenith angle. 

Finally the observed values of qQ which are used in con- 
nection with the solution of (42) for q^ should be carefully 
corrected for initial calibration and pyrheliometer degrada- 
tion errors after the procedure recommended by Hanson et al. 
(1972). 
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APPENDIX TABLE A 



NORTH AMERICAN PYRHELIOMETER STATION INFORMATION 



St&ti 



No. 


Code 


Name 


Index 


Lat . 


Long. 


Albedo 


1 


AA 


Fort Worth, Texas 


72259 


32-50 


97-03 


16 


2 


AB 


Apalachicola, Fla. 


220 


29-44 


84-59 


12 


3 


AC 


Miami, Florida 


202 


25-48 


80-16 


15 


4 


AD 


Brownsville, Texas 


250 


25-54 


97-26 


16 


5 


AE 


Charleston, S.C. 


208 


32-54 


80-02 


13 


6 


AF 


Nashville, Term. 


327 


36-07 


86-41 


14 


7 


AG 


Columbia, Mo. 


445 


38-58 


92-22 


16 


8 


AH 


Dodge City, Kansas 


451 


37-46 


99-58 


18 


9 


AI 


Caribou, Maine 


712 


46-52 


68-01 


14 


10 


AJ 


Madison, Wise. 


64l 


43-08 


89-20 


17 




RS 


Green Bay, Wls. 


645 








11 


AK 


El Paso, Texas 


270 


31-48 


106-24 


19 


12 


AL 


Albuquerque, N.M. 


365 


35-03 


106-37 


22 


13 


AM 


Ely, Nevada 


486 


39-17 


114-51 


19 


14 


AN 


Los Angeles, Ca. 


295 


33-56 


118-23 


15 




RS 


San Diego, Ca. 


290 








15 


AO 


Phoenix, Arlz. 


278 


33-26 


112-01 


22 




RS 


Tucson, Arlz. 


274 








16 


AP 


Santa Maria, Ca. 


394 


34-54 


120-27 


15 




RS 


Vandenburg, Ca. 


393 








17 


AQ 


Bismark, N.D. 


764 


46-46 


100-45 


15 


18 


AR 


Great Falls , Mont . 


775 


47-29 


111-22 


17 


19 


AS 


Medford, Oregon 


597 


42-22 


122-52 


16 


20 


AT 


Tacoma, Wash. 




4 7-57 


124-33 


16 




RS 


Tatoosh, Wash. 


798 








21 


AU 


Fresno, Ca. 


389 


36-46 


119-43 


13 




RS 


Oakland , Ca . 


493 








22 


AV 


Cape Hat t eras, N.C. 


304 


35-16 


75-33 


14 


23 


AW 


Sterling, Va. 




39-00 


77-26 


15 


24 


AX 


Boston, Mass. 


509 


42-22 


71-02 


15 




RS 


Nantucket, Mass. 


506 








25 


AY 


New York, N.Y. 


74486 


40-39 


73-47 


15 


26 


AZ 


Omaha, Nebr. 


72553 


41-22 


96-03 


15 
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Stat . 



No. 


Code 


Name 


Index 


Lat . 


Lone . A 


lbedo 


27 


CA 


Beaverlodge CDA, Alta. 




55-12 


119-30 


14 




RS 


Prince George B.C. 


72896 








28 


CB 


Churchill, Man. 


913 


58-45 


90-05 


15 


29 


CC 


Departure Bay, B.C. 




49-12 


124-00 


15 




RS 


Tatoosh, Wash. 


798 








30 


CD 


Edmonton, Alta. 


74119 


53-19 


113-35 


15 


31 


CE 


Fredericton CDA, N.B. 




45-52 


66-32 


15 




RS 


Caribou, ME 


72712 








32 


CF 


Goose A, (Nfld) 


816 


53-19 


60-25 


16 


33 


CG 


Guelph OAC (Ont.) 




43-36 


80-18 


14 




RS 


Buffalo, N.Y. 


528 








34 


CH 


Halifax, N.S. 


600 


44-53 


63-31 


16 




RS 


Sable, ISC 








35 


Cl 


Kapuskasing A. Ont. 


836 


49-25 


82-28 


15 




RS 


Koosonee CDA, Ont. 








36 


CJ 


Kentville CDA N.S. 




45-06 


64-30 


15 




RS 


Portland, ME 


606 








37 


CK 


Montreal, Qua. 




45-30 


73-42 


16 




RS 


Maniwaki, CDA, Que 


722 








38 


CL 


Normanain, CDA, Que 




48-54 


72-36 


15 




RS 


Sept. Isles, CDA, Que 


811 








39 


CM 


Ottawa NRC , Ont . 




45-24 ' 


75-42 


16 




RS 


Maniwaki , CDA, Que 


722 








4o 


CN 


St. John's West , CDA , (Nfld ) 




47-37 


52-45 


16 


4l 


RS 


Stephenville , (Nfld) 


815 


50-16 






CO 


Suf field A, Alta. 




111-11 


15 




RS 


Great Falls, Mont. 


775 








42 


CP 


Summer land CDA, B.C. 




49-36 


119-42 


16 


43 


RS 


Spokane, Wash. 


785 




107-41 




CQ 


Swift Current CDA, Sask. 


768 


50-17 


15 


44 


RS 


Glasgow, Mont. 


43-43 






CR 


Toronto , Ont . 




79-14 


14 




RS 


Buffalo, N.Y. 


528 








45 


CS 


Toronto Scarborough, Ont. 




43-43 


79-14 


14 




RS 


Buffalo, N.Y. 


528 








46 


CT 


Vancover U3C, B.C. 




49-11 


123-10 


16 




RS 


Tatoosh, Wash. 


798 








47 


CU 


Winnipeg Int A. Man. 


7 47 


49-54 


97-14 


13 




RS 


International Falls, Minn. 
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APPENDIX B 



Tables of observational data Interpolated to the effec- 
tive scan-time at each of the 47 stations used In this study, 
and for days 1, 5* (July 15 19. respectively). 

The symbols Q s , Q R , Qq , Q A , are as defined In Eq . 15 
and In the paragraph which Immediately follows this equa- 
tion. The symbol U denotes the preclpltable water vapor 
at the station and CZ Is the cosine of the zenith angle at 
the effective time of the observation. 

DELETE CODE symbols P and A Indicate Inconsistent trans- 
mission and absorption values, respectively, and these data- 
cases were not used. All other codes, such as A MAX, A HIM, 

T MIN express the occurrence of statistically Improbable 
values, which cases were used In a secondary data-cases 
(Sec. III. (C)). 
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